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LECTURE 22. CONVOLUTION 



Motivation: buildup of a pollutant in a lake. Let's say we have a lake and a pollutant is being 
dumped into it at the variable rate f{t). The pollutant degrades over time exponentially. If the 
lake begins at t = with no pollutant, how much is in the lake at time t > 0? 

The small drip of pollutant added to the lake between ti and ti+At, where At small, is f{ti)At. 
Later t > ti, the drip reduces to e~^^^^^^^ f{ti)At, where a > is the decay constant. Adding them 
up, starting at the initial time ti = 0, we obtain that the amount is 

(22.1) /*e-"(*-*i)/(ti)dti. 



Integral of this kind is called a convolution. 

We can solve this problem by setting up a differential equation. Let y{t) be the amount of 
pollutant in the lake at time t. Then, the amount of the chemical in the lake at time t + At is the 
amount at time t, minus the fraction that decayed plus the amount newly added: 

y{t + At) = y(t) - ay(t)At + /(t)At. 

Taking the limit as At ^ we obtain 

y' + ay = f{t), y(0) = 0. 

It is straightforward that (22.1) gives the solution of the above initial value problem. 

The convolution integral. The convolution of / and g is defined as 

(22.2) {f*9m= [' f{h)g{t-h)dh. 

Jo 

It gives the response at the present time t as a weighted superposition over the input at times 
ti < t. The weight g{t — ti) characterizes the system and /(ti) characterizes the history of the 
input. To ensure the existence of the integral, in what follows, we assume that f,g £ A. 

Example 22.1. Let /(t) = e^^* and g{t) = e^^*, where Bi ^ B2 are constants. Then, 



if*9m 



Jo 



Bi — B2 

Since £e"* = l/(s — a) for any constant a, we have 

It is not a coincidence, but rather, it is a property of convolution, as discussed below. 

The convolution operator acts like ordinary multiplication in that, if /, g, h are admissible then 

(i) (distibutive) f*{g + h) = f*g + f*h, 

(ii) (commutative) f * g = g * f, 

(iii) (associative) f * {g * h) = {f * g) * h. 
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However, the convolution operator differs from the multipHcation operator. For example, f*l ^ f, 
f*f^fiT^ general. 

Exercise. Show that * 1 = and cos t * cos t = ^{t cos t + sin t). 

Nevertheless, convolution in the t-domain does corresponds to multiplication in the s-domain. 

Theorem 22.2 (Convolution Theorem). If f,g e E, then f * g e E and C{f * g) = {£f){Cg). 

Proof. If f,g £ E, then f * g is continuous for all t £ [0,oo). Since \ f{t)\ ^ Aie^^* and \g{t)\ ^ 
^26^2*, then 



\{f*9m\ ^ rAie^^%e^^(*'*^)dti = ^1^2^^ 
Jo ^1 



Bit _ ^B2t 



B2 

Therefore, f * g £ E. Better yet, |/| * l^l € E. In other words, £(/ * g) converges absolutely for 
large s. 

For simplicity, let f{t) = and g{t) = for t < 0, so that 

/OO f'OO fOO 

e-''fit)dt, Cg= e~''g{t)dt, and {f * g){t) = f{t^)g{t-t^)dh. 
-OO J OO J~oo 

Indeed, the lower limit of the integral of f * g can be replaced by — oo since /(ti) = for ti < 
and the upper limit can be replaced by oo since g{t — ti) = for t — ti < 0. Then, 



/OO /"OO 
f{ti)g{t-ti)dhdt 
-oo J ~oo 

/oo / I'OO \ 
/ e-^*g(t - ti) dt f{h) dti 
-oo \J —oo / 



oo / fOO \ 

e-^(*^+*2)(7(i2) dt2 fih) dh 

— oo — oo / 

/oo 
e-''^Cg{s)f{h)dh = {Cf){Cg). 
-oo 



□ 



Convolution allows an easy passage from the s-domain to the t-domain and leads to explicit 
solutions for a general inhomogeneous term f{t). 

Example 22.3. Solve the initial value problem 

y" + co^y = comfit), y(0) = y'(0) = 0, 
where is constant and f £ E. 
Solution. Taking the transform. 

Since C( sin ujt) = -^r^ — ^, it follows that Ly = (£/)(£cjsinu;t). By the convolution theorem 
and by uniqueness, then 

= / * (w sinoit) = a; / f{t\)'svQ.uj{t — t\)dt\. 
Jo 

Note that we have a formula for the rest solution corresponding to the arbitrary function /. 
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The tautochrone. Starting from the rest state, suppose that a particle of mass m slides down a 
frictionless curve under gravity as shown in the figure below. 

The aim is to determine the shape of the curve so that the 
time of descent is independent of the starting point. A curve 
of this kind is called a tautochrone. It comes from two Greek 
words, meaning "same" and "time." The problem was solved 
by the Dutch mathematician Christian Huygens in 1673 as part 
of his theory of pendulum clocks. 

Let y be the height at which the particle starts and let v be the 
speed of the particle at the height z. The change in the kinetic 
energy (l/2)mi;^ at the height z is the change in the potential 
energy, which leads to 




Figure 22.1. The tautochrone. 



(22.3) 
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where g is the gravitational acceleration. Let a = a{y) be the arc from the rest state to the lowest 
point. The time of descent is 



y 1 da 

V dz 



dz 



y 1 
V 



(z) dz. 



da/dy\ 



y=z- 



Since the time of decent is constant, in view of 



Here, (j){y) = da/dy so that (j){z) 
(22.3), the problem reduces to 

[%{z){y - z)-^l^ dz = ci, 
Jo 

where ci is constant. The left side is the convolution of (j) and y^^l"^ . Taking the transform, then 



y-i/2) = (£,^)(£y-V2)=£ci 



ci 
s 



The first equality uses the convolution theorem. 

Since* li\t^^l'^\ = ^/vrs"^/^, we have 

-1/2 



-1/2 



£0 = C2S 0(y) = cy 

where c, C2 are constants. Since = da/dy, the equation of the curve reduces to 



c 

y 



or dx 



Idy. 



Exercise. By parameterizing y = sin^ 9, show that the curve is a cycloid 



-(26* + sin 26*), 



^1 - cos 2(9). 



When a particle slides down a frictionless curve, a related question is to find a path which 
minimizes the time of descent. The minimizing curve is called a "brachistochrone". It is interesting 
that the solution is again a cycloid, (not a straight line!) 



See the pset question on the Gamma function. 
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